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ABSTRACT 

Supersonic flows are expected to exist in the atmospheres of irradiated exoplanets, but the question of whether 
shocks develop lingers. Specifically, it reduces to whether continuous flow in a closed loop may become su- 
personic and if some portions of the supersonic flow steepen into shocks. We first demonstrate that continuous, 
supersonic flow may exist in two flavors: isentropic and non-isentropic, with shocks being included in the 
latter class of solutions. Supersonic flow is a necessary but insufficient condition for shocks to develop. The 
development of a shock requires the characteristics of neighboring points in a flow to intersect. We demon- 
strate that the intersection of characteristics may be quantified via knowledge of the Mach number. Finally, we 
examine 3D simulations of hot Jovian atmospheres and demonstrate that shock formation is expected to occur 
mostly on the dayside hemisphere, upstream of the substellar point, because the enhanced temperatures near 
the substellar point provide a natural pressure barrier for the returning flow. Understanding the role of shocks 
in irradiated exoplanetary atmospheres is relevant to correctly modeling observables such as the peak offsets of 
infrared phase curves. 

Subject headings: planets and satellites: atmospheres 



1. INTRODUCTION 

Current astronomical techniques favor the detection of 
close-in, irradiated exoplanets that are believed to be tidally 
locked. With temperatures reaching ~ 1000-3000 K, the 
flows in these atmospheres are expected to approach or exceed 
the local sound speed (~ 1 km s _1 ). A lingering question in 
the astrophysical literature is whether these supersonic flows 
can develop into shocks and subsequently dissipate some of 
the kinetic energy of the flow as heat. Specifically, the ques- 
tion reduces to whether continuous flow in a closed loop may 
become supersonic and if some portions of the supersonic 
flow steepen into shocks. 

In the present Letter, we provide the answers to these ques- 
tions by re-visiting a problem that is mathematically and 
physically similar: the design of wind tunnels and the de 
Laval nozzle. In $2] we demonstrate the possibility of super- 
sonic flow in a closed-loop flow and state a simple criterion 
for determining if these flows develop shocks. In S|3] we ap- 
ply our method to simulations of hot Jovian atmospheres. The 
implications of our study are discussed in S]4] 

2. THE DE LAVAL NOZZLE AND WIND TUNNEL DESIGN 

While much of thi s sec tion is based on the monograph of 
Liep mann & Roshkol (119571) . the re-derivation of these results 
elucidates the key concepts and assumptions made in stating 
our method, which we will later apply to hot Jovian atmo- 
spheres. 

2.1. Demonstrating the Possibility of Supersonic Flow 

In its most basic form, the de Laval nozzle consists essen- 
tially of a tube that continuously narrows to a minimum in its 
cross section, known as the "throat", before widening again in 
a symmetric fashion. The conservation of mass dictates that 
the mass flux M = pvA is constant, such that 

dA dv ^ dp 
A v p 
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where A is the cross-sectional area of the nozzle, v is the flow 
speed normal to the cross section and p is the mass density 
of the fluid. Invoking the conservation of momentum (dP = 
pvdv with P being the pressure) and the definition for the 
sound speed (c? s = dP/dp), we obtain 



dp = _ M2 d» t 
p v 



(2) 



with Ai = v /c s being the Mach number of the flow. When 
the flow is strictly subsonic (Ai <C 1), it becomes virtually 
incompressible (dp « 0). Substituting this expression into 
the equation for the conservation of mass yields 



dA dv - 2 , 
A v v ; 



(3) 



When the flow is subsonic (Ai < 1), narrowing the nozzle 
(dA < 0) leads to an increase in the flow speed (dv > 0). The 
flow is symmetric with respect to the throat. However, if the 
flow is supersonic (Ai > 1), narrowing the nozzle leads to the 
somewhat counter-intuitive result that a decrease in the flow 
speed (dv < 0) occurs. The flow is now asymmetric across 
the throat because of compressibility effects. In the design of 
continuous wind tunnels there is often a region of narrowed 
cross section, downstream of the de Laval nozzle, known as a 
"diffuser", which decelerates the flow and discharges it back 
into the return circuit. When dA = 0, we get dv ^ only 
when M = 1. 

The entry into the de Laval nozzle is described by a (small) 
velocity v m , a temperature T ln , a pressure P m and a mass den- 
sity p m . By enforcing the conservation of energy and assum- 
ing an ideal gas, we may transform the isentropic relations 
into ones for the pressure and density ratios, 
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Furthermore, the Mach number is reduced, 



A s /A 

FIG. 1 . — Subsonic and supersonic solutions for the de Laval nozzle as func- 
tions of A s I A, where A is the cross-sectional area and A s is the cross- 
sectional area corresponding to the sound speed. Top panel: Mach num- 
ber. Bottom panel: ratio of pressure to entry pressure. The sonic point is 
marked by a filled circle. The empty/unfilled circles indicate a possible, non- 
isentropic solution for a M = 2 shock. 



where 7 = 1 + 2/ndof = 7/5 is the adiabatic gas index and 
fMof = 5 is the number of degrees of freedom in the gas. Note 
that Pj n > P and p m > p. By enforcing the conservation of 
mass, one obtains the cross-sectional area normalized by its 
sonic value, 
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By utilizing equation ©, one can eliminate M. in favor of 

P/Plu, 
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Figure[T]shows M. and P/p n as functions of A s /A for 7 = 
7/5. The cross-sectional area A = A(x) serves as a proxy for 
the actual distance x along the nozzle. There are generally two 
solution branches: subsonic and supersonic. If a given nozzle 
has a throat with a cross-sectional area larger than A s , then 
the solution remains in the subsonic branch and is symmetric 
about the throat. If the throat has a cross-sectional area A s , 
then the system develops an asymmetry, where the flow from 
the entry point to the throat is subsonic, while the flow from 
the throat to the exit point is supersonic. The pressure drops 
monotonically from a value of p n to P ou t- If the flow remains 
isentropic, then the exit pressure P out is determined by the 
supersonic solution to the PJ P m curve in FigureQ] Generally, 
isentropic solutions cannot exist between the two bracketing 
values of P/p n for a given A s j A. However, non-isentropic 
solutions may exist. In the example shown, a Mach-2 shock 
occurs in the flow. The existence of the shock creates entropy 
and re-compresses the flow, 
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Post-shock quantities are indicated with a prime. 

While our present study of the de Laval nozzle cannot be di- 
rectly applied to exoplanetary atmospheres and thus does not 
yield any quantitative predictions, we do gain some heuristic 
insight: 

• We expect supersonic flow to be able to develop in 
continuous flows inherent in exoplanetary atmospheres 
(without orography), consis tent with the results of pub- 
lished 3D simulations (see c]4. ll for references); 

• In the de Laval nozzle, variations in the pressure, 
density and Mach number are induced by changes in 
the cross-sectional area. In irradiated atmospheres, 
changes in the Mach number are induced by changes 
in temperature; 

• For close-in exoplanets that are expected to be tidally 
locked, the hemispheric (dayside-only) nature of the 
stellar irradiation naturally introduces re-compression 
into the continuous flow from day side to nightside. Irra- 
diated, tidally-locked exoplanets create their own pres- 
sure "barrier" for the atmospheric flow to crash into, 
thus creating shocks. We expect this pressure barrier to 
exist near the substellar point. 

2.2. Demonstrating the Existence of Shocks 

The question of whether a shock develops in a continuous 
flow can be reduced to analyzing how the Mach number of 
the flow changes along a given streamline. Since mass, mo- 
mentum and energy are conserved, the problem is mathemat- 
ically identical to employing the established machinery of the 
Rankine-Hugoniot jump conditions. 

Consider a flow in 2D, Cartesian coordinates. Let the flow 
be described by the velocity components v x and v y , such that 
v 2 = v 2 +v 2 . Across an interface, the Mach number changes. 
Let the angle between the flow and the interface be (3 such 
that tan/3 = v x /v y . Downstream of the interface, the flow 
velocity is v' and the flow direction is deflected by an angle 
8 relative to v, yielding tan(/3 — 6) = v' x /v'. Denoting the 
upstream Mach number by A4 = v/c s , the ratio of the down- 
stream to the upstream mass density is 



£ _ (7 + 1) M 2 sin 2 p 
~~p ~ (7- 1) M 2 sin 2 13 + 2' 



(9) 



A relationship between (3 and 9 may be obtained by taking the 
ratio of tan j3 to tan(/3 — 9) and invoking the jump condition 
for p'/p, 



M 2 sin 2 f3 - 1 



(7 + 1) A^ 2 sin/3sin( 
2 cos {13-9) 



(10) 



For 6 = 0°, two corresponding values of f3 exist. The first is 
(3 = 90°, which describes a situation where the downstream 
flow is subsonic (e.g., for a strong shock). The second occurs 
when 

A, = sin" 1 (1/M) , (11) 

which describes a situation where the downstream flow is su- 
personic. 
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FIG. 2. — Characteristic angles /3o, at the equator, computed using the 
Showman & Polvani (2011) analytical model. Shown are models with dif- 
ferent values of the dimensionless radiative (t ra( j) and hydrodynamic (tdrag) 
drag time scales. For the model with t ra( j = tdrag = 1. the dimensionless 
velocity (see text for definition) is always less than unity and thus /3o is un- 
defined. The dimensionless distance x is given in units of the Rossby length 
scale. 

Supersonic flow is a necessary but insufficient condition 
for a shock to develop. A wave with a constant propagation 
speed possesses the property that every point on the wave is 
translated by the same distance after a given time interval. 
In the parameter space of distance versus time, a given point 
on the wave at different times is connected by a straight line 
known as a characteristic. It is apparent that these character- 
istics never intersect. In reality, the propagation speed also 
depends on the density of the fluid, implying that different 
points on the wave are translated by different distances over 
a given time interval. In this case, the characteristics may 
intersect — if they do, a shock forms and converts part of the 
kinetic energy into heat. 

At the equator of an irradiated, tidally-locked exoplanet, 
the flow is predominantly zonal (east-west). For purely zonal 
flow, each characteristic is associated with its own value of 
the angle (3q (sometimes known as the "Mach angle"). When 
/?o decreases, the characteristics diverge. When (3q increases, 
they converge. Convergence is not necessary an indication 
of intersection — there needs to be enough room for the char- 
acteristics to intersect. On a sphere (with a radius R) where 
two points are separated by a small angular distance A<f>, the 
horizontal distance h over which intersection occurs is given 
by 



R \tan/?o,i tan/?o,2/ ' 

where the subscript "2" refers to a point downstream of an- 
other point subscripted by "1". 

The method stated in this sub-section may be applied to 
local pairs of points in a global flow. 

3. APPLICATION TO HOT JUPITERS 

3.1. Analytical Models 

We first apply the method described in the previous section 
to the 2D, linear, analyt i cal, 1 .5-layer, shallow water model 
of iShowman & Polvanil ( 1201 11) . which reproduces some of 
the flow features of hot Jovian photospheres. There are sev- 
eral important limitations associated with the shallow water 
model. The governing equation for the shallow water height 
derives from the condition of incompressibility. There is no 



governing equation for temperature. Both features imply that 
a sound speed cannot formally be defined. However, a proxy 
for the sound speed is the gravity wave s peed , which is de- 
fined in shallow water systems as c g = y 'gH where g is the 
surface gravity and H is the equilibrium water height. When 
the flow velocity exceeds c g , the analogues of shocks known 
as hydraulic jumps develop. The flow velocity normalized by 
c g is then a proxy for the Mach number. 

Stellar irradiation and radiative transfer are not explicitly 
considered in the shallow water model. Rather, the effects of 
radiation are mimicked using a dimensionless radiative drag 
time scale i ra( j. A dimensionless time scale idrag mimics 
the effect of hydrodynamic drag. Figure [2] shows calcula- 
tions corresponding to the 9 models presented in Figure 3 of 
IShowman & Polvanil d201 lh . We have shown the characteris- 
tic angle (3 = (3q only at the equator, since this simplifies our 
analysis. Negative values of f3o correspond to counter-rotating 
(rather than super-rotating) flow. 

First, we see that there are locations where the magnitude 
of /3q increases and thus hydraulic jump formation may oc- 
cur. Second, the existence of these locations is robust to the 
strength of irradiation (through < ra d) and the presence of drag 
(through idrag)- However, since we do not expect the shallow 
water model to completely reproduce the flow field computed 
from 3D simulations, we do not expect it to correctly mimic 
the locations of shock formation. 

3.2. 3D Simulations 

We next apply our method to 3D simulations of at- 
mospheric circulation of hot Jupiters, which solve for 
atmospheric dynamics and dual-band radiative trans- 
fer self-consistentlv faeng. Menou & Philhppsl 1201 It 
iHeng. Frierson & Phillippsll201 lh. S pecific ally, we first adopt 
Model H of iPerna. Heng & Ponl (120121) as a case study, 
which has an equilibrium temperature of T cq w 1720 K 
and has no temperature inversion present in its atmosphere 
(70 = 0.5). (See Tables 1 and 2 of lPerna. Heng & Pontll2012l 
for more details.fl 

The top panel of Figure [3] shows the temperature-velocity 
map near the infrared photosphere of our model hot Jupiter. 
As witnessed in published 3D simulations of hot Jovian at- 
mospheres, a global, chevron-shaped feature is centered about 
the equator, which is a mostly linear response arising from the 
interac tion of the mean flow with s tanding Rossby and Kelvin 
waves (Showm an & Polvanill201 ll) . 

The middle panel of Figure[3]shows the characteristic angle 
/3o versus longitude at the equator. The corresponding Mach 
numbers are supersonic throughout. In the direction of the 
zonal flow, flo increases with longitude largely on the dayside 
hemisphere. At locations where shock formation has a chance 
of occurring, we estimate the values of h/R (bottom panel 
of Figure [3$. In the case of Model H, locations exist where 
h/R < ir/2 implying that the intersection of characteristics, 
associated with points in the flow where (3q is increasing, will 
occur. We reach the same conclusion when we examine atmo- 
spheric layers at higher altitudes (lower pressures). At pres- 
sures of P > 1 bar, the equatorial Mach numbers become less 
than unity and the necessary condition for shock formation is 
unfulfilled. 

In F igure 2J we examine Model W of lPerna. Heng & Pontl 
(120121) . which has T cq rj 970 K. Our conclusions drawn 

3 We note a typographical error in Table 1 of Perna, Heng & Pont 12012): 
rg should have values of 5 X 10 3 and 2 X 10 4 . 
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FIG. 3. — Determining the locations where shocks form near the infrared 
photospheres of hot Jupiters. The equilibrium temperature adopted in this 
model is about 1720 K (Model H). Top panel: map of temperature (colors; 
in units of K) and velocity (arrows) at P f» 10 mbar. Middle and bottom 
panels: /3rj and h/R, at the equator, as functions of longitude. 

from examining Model H carry over: shocks form mostly 
on the dayside hemisphere, where fio is an increasing func- 
tion of longitude. We also examine Models W and H with 
temperature inversions present (70 = 2). Besides minor 
quantitative differences between the results from these mod- 
els and those presented in Figures[3]and[4] we find no qualita- 
tive differences. Models with lower equilibrium temperatures 
(T eq w 540-750 K; Models C, CI and C2) do not develop 
supersonic flows. 

In all of the cases examined, shocks — when they do form — 
are expected to form mainly upstream of the substellar point. 
The enhanced temperatures near the substellar point are cre- 
ating a natural pressure barrier for the returning flow, from 
the nightside, to "crash " into, thus creating conditions con- 
ducive to shock formation^ While we have focused on hot 
Jupiters, we expect the insight gained from the present study 
to apply to all types of hot exoplanets — Earth-like, Neptune- 
like or Jupiter-like. The atmospheres of tidally-locked, irradi- 
ated exoplanets are akin to forced wind tunnels. 

4. DISCUSSION 

4 We may gain further insight by recalling the classical situation of a 
M 3> 1 flow impinging upon a physical barrier, which forces the down- 
stream velocity to ultimately become zero. In this case, an increasing 0q is 
enforced by decreasing v. In irradiated exoplanetary atmospheres, an increas- 
ing /3q is enforced by increasing T. 
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FIG. 4. — Same as Figure|3] but for an equilibrium temperature of about 970 
K (Model W). 

4.1. Relevance to Previous Work 

Previously published work on simulations of irradiated ex- 
oplanetary atmospheres have recorded supersonic flows when 
the temperature reaches ~ 1000 K. Most of these studies 
utilize "general circulation models" (GCMs), which are tradi- 
tionally designed to model the terrestrial cli mate system and 
do no t include an explicit treatmen t of shocks ( Showma n et al. 
| 2009t iThrastarson & Choi 1201 U iHeng Menou & Phillipps 
201 it IHeng. Frierson & Phillipps! 1201 U iRauscher & Menou 
120121) . These simulations may describe supersonic flow, but 
do not correctly convert a fracti on of the kinetic energy int o 
heat. Exception s a re the wor k oflDobbs-Dixon & Lin! (1 2008). 
iLi & Goodman! (120101). [D obbs-Dixon, Cumming & Lin 
d2010l) and iDobbs-Dixon. Agol & Burrows! d20 121) . which 
include proper shock capturing. The simulations of 
Li & Goodman (2010) are 2D and require confirmation by 
3D simulations due to the global nature of flow features in 
the atmospheres of close-in, irradiated e xoplanets. Although 
their simu lations are 3D in nature, Dobbs-D ixon & Lin! 
(2008) and [Dobbs-Dixon. Cum ming & Lin (2010) employed 
a non-glob al grid in an attempt to c ircumvent the "pole 
problem" (Staniforth & Thuburn 2012), where the conver- 
gence of meridians at the poles results in vanishing time 
steps for the simulations. Dobbs-Dixon, Agol & B urrows! 
d2012l) employed a partial solution to the pole problem that 
vio lates causality. Detailed compa r isons between the work 
of IDobbs-Dixon. Agol & Burrows! d2012l) and the insight 
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gained from the present study require future 3D simulations 
that solve the pole problem and include shock heating in a 
s elf-consistent, energy - conse rving manner. 

iPerna. Heng & Pong (120 121) examined the effects of varying 
stellar irradiation on the atmospheric circulation. At locations 
where the Mach number exceeds unity, they assumed a shock 
to form and dissipate a fraction 4(7 — l)/(7 + l) 2 of the ki- 
netic energy into heat. In light of our present results, such an 
approach produces an over-estimation of the shock heating. 
Nevertheless, it remains likely that shocks do not penetrate 
deeply enough into the atmosphere to affect the evolution of 
the exoplanet in a significant manner. 



4.2. Observational Consequences 

We have demonstrated that shocks are expected to be 
present near the infrared photospheres of sufficiently irradi- 
ated exoplanets. The main effect of shocks is to convert some 
fraction f s of the local kinetic energy of the atmosphere into 
heat. We now wish to derive an expression for f s and also es- 
timate the ratio r s of the flux of shock heating to the incident 
stellar flux. 

We start from the shock jump condition for temperature 
(Liepmann & Roshko 1957), 



2kT'/mv 2 and f So — 2kT/mv 2 , we obtain 



fs J so 



4(7-1) 
(7 + I) 2 



1 - 



1 



(7-1) 

7X 2 jM 4 



(14) 



where f Sa is the fraction of thermal to kinetic energy in the 
pre-shock flow. For example, if f So = and M. = 2 then 
we obtain f s w 0.25. The ratio of this shock heating to the 
incident stellar energy is 



fs<JsBk B v 2 T 3 
2Tq Kcpm 



exp(rs), 



(15) 



where ctsb denotes the Stefan-Boltzmann constant, fee the 
Boltzmann constant, t§ ~ 1 the shortwave optical depth, 
the top-of-the-atmosphere stellar flux, n — 2/(2 + ridof) the 
adiabatic coefficient, cp the specific heat capacity at constant 
pressure and m the mean molecular mass. For Model H, we 
estimate that r s ~ 0.1 at locations where shock formation 
occurs. 

Shock heating functions as an additional form of drag, 
which acts to reduce the efficiency of heat redistribution from 
the dayside to the nightside hemisphere of a exoplanet. Con- 
sequently, the day-night flux contrast is expected to increase, 
while the shift of the peak of the phase curve from the sub- 
stellar point is expected to decrease. Simulations that fail to 
capture shock heating will thus over-estimate the peak offset. 



Using 7.M 2 = mv 2 /kT and recognizing that f s 
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